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Method of Seismic Signal Processing and Exploration 



Technical Field 

This invention relates to the general subject of seismic exploration 
5 and. in particular, to methods for identifying structural and stratigraphic 
features in three dimensions. 



Background Of The Invention 
Ordinary 2-D seismic data is acquired along lines (See lines 10 and 

10 11 in FIG. 1) that consist of geophone arrays onshore or hydrophone 
streamer traverses offshore. Geophones and hydrophones act as sensors to 
receive energy that is transmitted into the ground and reflected back to the 
surface from subsurface rock interfaces 12. Energy is usually provided 
onshore by vibroseis vehicles which transmit pulses by shaking the ground 

15 at pre-determined intervals and frequencies on the surface. Offshore, airgun 
sources are usually used. Subtle changes in the energy returned to surface 
often reflect variations in the stratigraphic, structural and fluid contents of the 
reservoirs. 

In 3-D seismic the principle is similar, however, lines and arrays are 
20 more closely spaced (See FIG's. 1 and 2) to provide more detailed sub- 
surface coverage. With this high density coverage, extremely large volumes 
of digital data need to be recorded, stored and processed before final 
interpretation can be made. Processing requires extensive computer 
resources and complex software to enhance the signal received from the 
25 subsurface and to mute accompanying noise which masks the signal. 

Once the data is processed, geophysical staff compile and interpret 
the 3-D seismic information in the form of a 3-D cube (See FIG. 4) which 
effectively represents a display of subsurface features. Using the data cube, 
information can be displayed in various forms. Horizontal time slice maps 
30 can be made at selected depths (See FIG. 5). Using a computer workstation 
an interpreter can slice through the field to investigate reservoir issues at 
different horizons. Vertical slices or sections can also be made in any 
direction using seismic or well data. Time maps can be converted to depth 
to provide a structural interpretation at a specific level. 
35 Three-dimensional (3-D) seismic is being used extensively worldwide 

to provide a more detailed structural and stratigraphic image of subsurface 
reservoirs. Acceptance of 3-D seismic has accelerated during the last five 
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years based on a proven track record that continues to grow. The 3-D 
payout has been measured by increased reserve estimates, cost savings 
from more accurate positioning of delineation and development wells, 
improved reservoir characterization leading to better simulation models, and 
5 the ability to predict more accurately future opportunities and problems 
during the production history of a field. More importantly, 3-D seismic has 
also been used as an exploration tool to reduce drilling risk in structurally 
complex areas and to predict reservoir quality in undrilled areas. 

As good as 3-D seismic surveys and interpreters have become, 
1 0 improvements are needed. 

In particular, seismic data has been traditionally acquired and 
processed for the purpose of imaging seismic reflections. Changes in 
stratigraphy are often difficult to detect on traditional seismic displays due to 
the limited amount of information that stratigraphic features present in a 
15 cross-section view. Although such views provide an opportunity to see a 
much larger portion of these features, it is difficult to identify fault surfaces 
within a 3-D volume where no fault reflections have been recorded. More 
importantly, seismic data is not known to have been acquired or used for the 
purpose of imaging seismic discontinuities instead of seismic reflections. 

20 

Summary Qf The Invention 
In accordance with the present invention, a method is disclosed for 
the exploration of hydrocarbons. The method comprises the steps of: 
obtaining a set of seismic signal traces distributed over a pre-determined 

25 three-dimensional volume of the earth; dividing the three-dimensional 
volume into a plurality of vertically stacked and generally spaced apart 
horizontal slices; dividing each of the slices into a plurality of cells that are 
arranged into laterally extending rows and columns and that have portions 
of at least three generally vertically extending seismic traces located therein; 

30 measuring across each of the cells the cross-correlation between one pair of 
traces lying in one vertical plane to obtain an in-line value and measuring 
the cross-correlation between another pair of traces lying in another vertical 
plane to obtain a cross-line value that are estimates of the time dip in an in- 
line direction and in a cross-line direction; combining the in-line value and 

35 the cross-line value to obtain one coherency value for each of the cells; and 
displaying the coherency values of the cells across at least one horizontal 
slice. 
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This technique is particularly well suited for interpreting fault planes 
within a 3-D seismic volume and for detecting subtle stratigraphic features in 
3-D. This is because seismic traces cut by a fault line generally have a 
different seismic character than traces on either side of the fault. Measuring 
trace similarity, (i.e., coherence or 3-D continuity) along a time slice reveals 
lineaments of low coherence along these fault lines. Such coherency 
values can reveal critical subsurface details that are not readily apparent on 
traditional seismic sections. Also by calculating coherence along a series of 
time slices, these fault lineaments identify fault planes or surfaces. 

Numerous other advantages and features of the present invention will 
become readily apparent from the following detailed description of the 
invention, the embodiments described therein, from the claims, and from the 
accompanying drawings. 

Brief Description of the Drawing 
FIG. 1 illustrates an arrangement of geophones to obtain 3-D seismic 
data from the earth's subsurface for processing in accordance with the 
present invention; 

FIG. 2 is a plan view of the arrangement shown in FIG. 1 ; 
FIG. 3 is a representation of the seismic traces laying in a plane 
passing through one row of geophones shown in FIG. 2; 

FIG. 4 is a pictorial representation of the information obtained from 
processing 3-D seismic data; 

FIG. 5 is a pictorial representation of a horizontal time slice of 3-D 
25 seismic data processed in accordance with the prior art; and 

FIG. 6 is a pictorial representation of a horizontal time slice of 3-D 
seismic data processed in accordance with the present invention. 

Detailed Description 

30 While this invention is susceptible of embodiment in many different 

forms, there is shown in the drawings, and will herein be described in detail, 
one specific embodiment of the invention. It should be understood, 
however, that the present disclosure is to be considered an exemplification 
of the principles of the invention and is not intended to limit the invention to 

35 the specific embodiment or algorithm so described. 

The first step is to obtain a set of seismic data in the form of seismic 
signal traces distributed over a three dimensional volume of the earth. 



10 
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Methods by which such data is obtained and reduced to digital form for 
processing as 3-D seismic data are well known to those skilled in the art. 

The next step is to generate a "discontinuity cube." This is done by 
applying a coherency algorithm to the 3-D seismic data. This algorithm may 
5 take many forms. Whatever its form, its function is to compare the similarity 
of nearby regions of seismic data within the 3-D seismic volume. If a trace 
segment is similar to its neighbors (e.g., in the in-line and cross-line 
directions), it is assigned a low discontinuity value; if a trace segment is not 
similar to its neighbors, it is assigned a high discontinuity value. 
10 FIG. 2 is a plan view of a portion of 3-D seismic volume. In order to 

measure discontinuity, a trace segment at one point A is compared to 
adjacent trace segments B and C. One way to compute trace similarity is 
described below. 

The zero mean lagged cross-correlation in the in-line (x-direction) 
15 between trace u(t, x, y) and u(t, x+dx, y) with a lag time of "flag" msec is 
defined as: 

k = +w 

X u(t+k, x, y) u(t+k+tlag, x+dx, y) 
Px (t,tlag) = ^-^ 

Va(t, x, y) a(t, x+dx, y) 

20 where: 

k = +w 

a(t,x,y) = X u 2 (t+k, x, y) 
k = -w 

and 

k m +W 

25 a(t, x+dx, y) « X u 2 (t+k, x+dx, y) 

k * -w 

are autocorrelations used to normalize the cross-correlation, and where 
w+w is the length in msec of the correlation window. It is important to 
30 choose w large enough so the assumption of zero mean is valid. Values on 
the order of a seismic wavelet are appropriate. Other methods of 
normalization may be used (e.g., product of the energies of the traces, etc.). 
In particular, cross correlation is one method of combining two waveforms to 
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measure the similarities of the waveforms. Autocorrelation is a method of 
combining a waveform with itself. See Sheriff's "Encyclopedic Dictionary of 
Exploration Geophysics," Society of Exploration Geophysicists, Tulsa, 
Oklahoma. 

5 The zero mean lagged cross-correlation in the cross-line (y-direction) 

between trace u(t, x, y) and u(t. x, y+dy) with a lag time of tlag msec is 
defined as: 

k = +w 

X u(t+k, x, y) u(t+k+tlag, x, y+dy) 
P x (t, tlag) = 

Va(t, x, y) a(t, x, y+dy) 

where 

k = +w 

a(t, x, y+dy) = X u ( t+k « x - y +d y) 
k'= -w 

The direction of apparent time dip in the x and y directions is 
15 estimated to be that lag (i.e., tlagx and tlagy) that has the greatest (i.e., most 
positive) cross-correlation. These values are p x (t, tlagx) and p y (t, tlagy). 

Given the apparent dips (in msec/trace), it is a simple (but not 
necessarily accurate when dealing with noisy data) calculation to obtain dip 
and dip azimuth. More importantly, the concept of cross-correlation is 
20 extended to two dimensions by taking the geometric mean between the 
classical one dimensional cross-correlations: 

Pxy (t. tlagx, tlagy) = Vpx(t. tlagx) p y (t, tlagy) 

25 This value (or attribute) serves as a rather robust estimate of signal 

discontinuity within geologic formations as well as signal discontinuities 
across faults and erosional unconformities. 

Computer Pmpr m 

30 A simplified FORTRAN 77 program for performing these calculations 

is given below: 

Given a trace "x" from a 3-D seismic amplitude volume, and its two 
neighboring traces "y" (in the in-line direction) and "z" (in the cross-line 
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direction), subroutine COH calculates an output trace "rho M containing 
coherence coefficients using a running window cross-correlation algorithm 
where: 

"rnins" and w maxs K are the minimum and maximum sample indices for 
5 all four traces; 

"inwinr is the window length in samples; 

"nlags 1 ' specifies the number of lags (relative time shifts) to do each 
side of "CT in the cross-correlation; and 

"sr" is the sample interval in ms. 

At each sample, subroutine CROSS calculates a series of normalized cross- 
correlation coefficients, returning the largest coefficients for each direction in 
"rhol" and "rho2'\ The time shift at which the maximum coefficients occur is 
returned in "tshfl" and' M tshf2"; these times are not used. Subroutine COH is 
called repeatedly, once for every trace in the input seismic amplitude 
volume, to produce a new 3-D data volume or "coherency cube" containing 
coherence coefficients. 



subroutine coh (x, y, z, rho, mins, maxs, iwinl, nlags, sr) 
real x(mins:maxs), y(mins:maxs), z(mins:maxs) 
real rho(minsrmaxs) 

ihwin = iwinl/2 

do j = mins+ihwin, maxs-ihwin 
k = j - ihwin 

call cross (x(k), iwinl, y(k). iwinl, nlags, sr, rhol , tshfl) 
call cross (x(k) t iwinl, z(k), iwinl, nlags, sr, rho2, tshf2) 
rho(J) = sqrt (rho1*rho2) 

enddo 

return 
end 



subroutine cross (x, nx, y, ny, lags, sr, peak, tshift) 
real x(0:nx-1 }, y(0:ny-1 ), sr t peak, tshift 
parameter (maxlags=l 28) 

BNSDOCID: <WO 9618915A1_L> 



10 



15 



20 



25 



30 



35 



WO 96/18915 



PCTAJS95/13644 



10 



20 



real g(-maxlags:+maxlags) 
double precision xx.yy 

nlags = max(0, min(lags f maxlags)) 

tshift = 0.0 

peak = 0.0 

xx = 0.0 

yy = 0.0 

ks = 0 

do ix = 0, nx-1 
xx = x(ix)**2+xx 

enddo 



1 5 if (xx .eq. 0.0) return 

do iy = 0, ny-1 

yy = y(iy)**2 + yy 

enddo 



if (yy .eq. 0.0) return 



do is s -nlags, + nlags 
g(is) = 0.0 
25 . do it = 0, nx-1 

if (it-is .ge. 0) then 
if (it-is .le. ny-1) then 

g(is) = g(is) + x(it)*y(it-is) 
endif 

30 endif 

enddo 

if(abs(peak) .It. abs(g(is))) then 
peak = g(is) 
kssts 

35 endif 
enddo 

tshift = ks*sr 
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peak = peak/sqrt (xx*yy) 

return 

end 



Landmark and GeoQuest interpretive workstations, for example, can 
be used to view and interpret faults and stratigraphic features by loading the 

10 discontinuity cube as a seismic volume. Visualization software (e.g., 
Landmarks^ SeisCube software) may be employed to rapidly slice through 
the discontinuity volume to aid in understanding complex fault relationships. 
Discontinuity displays can reduce interpretation cycle time when used in 
selecting which seismic lines to interpret, enabling the interpreter to work 

1 5 around faults and poor data areas. In addition, subtle stratigraphic features 
and complex faulting which are not readily apparent on traditional seismic 
displays can be rapidly identified and interpreted. FIG's. 5 and 6 are side by 
side comparisons of the same seismic information displayed and processed 
conventionally and in accordance with the present invention. Fault lines are 

20 readily apparent in FIG. 6. 

Coherency maps have been run on several 3-D surveys. At depths of 
reasonable data quality, approximately 90% of the faults can be readily 
identified. Faults were identified on coherency maps which were very subtle 
on seismic sections, but clearly present on the coherency maps because of 

25 the robustness of the method and the map perspective of fault patterns. 
Since coherency maps can be run on uninterpreted time slices, the present 
invention offers a means to greatly accelerate mapping of the structural 
framework and to reveal details of fault relationships which would otherwise 
be interpreted only through tedious fault picking. 

30 

Specific Examples 

2-D seismic coherence maps were generated along 
picked horizons and clearly identified shale diapirs in offshore 
Nigeria. 

35 In offshore Gulf of Mexico, the technique readily 

identified diapiric structures. 

On several coherence time slices, remarkable detail of 
stratigraphic features, such as abandoned river channels, mud 
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flows, and submarine canyons, was displayed. On seismic 
sections, these features were sometimes apparent but, in some 
cases, were unidentifiable even with close scrutiny. 
This is the first known method of revealing fault planes within a 3-D 
5 volume where no fault reflections have been recorded. Faults are often 
critical to the accumulation of oil. A fault may form a seal by cutting off a 
structural or stratigraphic feature so the oil is trapped against the fault. On 
the other hand, if the fault plane contains rubble that has not been 
cemented, it may form a conduit for fluids. This may allow the hydrocarbons 
10 to drift up the fault plane into the feature and be trapped in it or to escape 
from the feature by drifting up the fault plane out of it. 

Thus, fault lines can predict flow patterns in a reservoir and 
communication between injector and producing wells, for example. Seismic 
discontinuities can also provide the needed link to enable reservoir 
15 prediction between the wells and establish reservoir continuity and flow 
patterns across a field. 

Coherency mapping of 3-D seismic is an extremely powerful and 
efficient tool for mapping both structure and stratigraphy. The new method is 
particularly sensitive to any lateral variation in wavelet character and 
20 therefore is particularly sensitive to the common causes of lateral variations 
in the wavelet (i.e., fault displacement or stratigraphic variations). This 3-D 
method analyzes a time-slice or horizon based interval and measures the 
maximum of the normalized cross-correlation in the in-line and cross-line 
directions. 

25 

Further Analysis 

The discontinuity cube will clearly highlight fault planes as zones of 
high discontinuity. However, these zones may not be clearly in areas of 
lower signal-to-noise ratio. A method of enhancing these fault zones 

30 involves the application of a "median planar operator." Faults in the earth's 
subsurface generally express themselves as planes or surfaces. In the case 
of a curved fault surface, a series of small flat planes may be used to 
approximate the fault surface. In accordance with this aspect of the 
invention, a small planar operator is used to enhance (i.e., a "filter") the 

35 identification of subtle stratigraphic features. First, a small region of seismic 
data is selected around a center value. This region may be formed from a 
plurality of the cells used to form the "coherency cube." A small fault plane is 
then mathematically inserted into the region, and the median value of the 
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points within the plane is calculated for the dip and azimuth that best aligns 
with the zone of high discontinuity. This median value is then assigned to 
the center value of a new array. Next, the region of data is shifted (e.g., by 
one row) and the process is repeated until each point in the previously 
5 determined discontinuity cube has been analyzed as a center value. The 
end result is a completely new discontinuity cube with fault planes enhanced 
and noise and stratigraphic features (i.e., non-planar features) attenuated. 
These stratigraphic features may be separated by subtracting the new 
discontinuity cube from the old discontinuity cube without the planar filter 
10 application. 

From the foregoing description, it will be observed that numerous 
variations, alternatives and modifications will be apparent to those skilled in 
the art. Accordingly, this description is to be construed as illustrative only 

15 and is for the purpose of teaching those skilled in the art the manner of 
carrying out the invention. Other algorithms may be used to measure the 
similarity of nearly regions of seismic data or to generate the "discontinuity 
cube." Moreover, equivalent computations may be substituted for those 
illustrated and described. Also certain features of the invention may be used 

20 independently of other features of the invention. 

For example, stratigraphic features have been generally identified on 
time slices where dips were low; and consequently, the time window 
captured a narrow statigraphic section. In areas of higher dip, the method 
should work on picked horizons. Therefore, as a stratigraphic mapping tool, 

25 there is good reason to believe that new levels of detail can be mapped than 
previously, although this may require mapping of the horizon of interest. 

As another example, while coherence slice maps by themselves are 
very powerful mapping tools, when used in conjunction with reconnaissance 
mapping of amplitudes and dip maps, there is promise of a technological 

30 milestone in mapping effectiveness for the Gulf of Mexico or similar basins 
with 3-D seismic. It is believed that detailed mapping of structure and 
stratigraphy will be accelerated by mapping in a map view and less by 
traditional line by line picking. Interpretation in a map view of 
-reconnaissance" data offers significant improvement in quality and quantity 

35 of interpretation. 

Thus it will be appreciated that various modifications, alternatives, 
variations, and changes may be made without departing from the spirit and 
scope of the invention as defined in the appended claims. It is, of course, 
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intended to cover by the appended claims alt such modifications involved 
within the scope of the claims. 
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CLAIMS 

We claim: 

1. A method for the exploration of hydrocarbons, comprising the 
5 steps of: 

a) obtaining a set of seismic signal traces distributed over a 
pre-determined three-dimensional volume of the earth; 

b) dividing said three-dimensional volume into a plurality of 
vertically stacked and generally spaced apart horizontal slices and dividing 

10 at least one of said slices into a plurality of cells that are arranged into 
laterally extending rows and columns, each of said cells having portions of 
at least three seismic traces located therein, each of said portions of said 
traces generally extending vertically through said cells, and a first trace and 
a second trace in said cell lying in one generally vertical plane and a third 

15 trace and said first trace in said cell lying in another generally vertical plane 
that is generally at right angles to said one vertical plane; 

c) measuring across each of said cells the cross- 
correlation between said traces lying in said one vertical plane to obtain an 
in-line value and the cross-correlation between said traces lying in said 

20 another vertical plane to obtain a cross-line value that are estimates of the 
time dip in an in-line direction and in a cross-line direction; 

d) combining said in-line value and said cross-line value to 
obtain one coherency value for each of said cells; and 

e) displaying said coherency values of said cells across at 
25 least one of said horizontal slices. 

2. The method of Claim 1, wherein step (b) is performed on each 
of said horizontal slices; and where in step (e) said coherency values are 
displayed over successive horizontal slices. 

30 

3. The method of Claim 1, wherein step (c) includes the step of 
normalizing each in-line value and each cross-line value. 

4. The method of Claim 3. wherein said step of normalizing 
35 includes the step of obtaining the product of the energies of each pair of 

traces. 
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5. The method of Claim 3, wherein the auto-correlation of said 
traces lying in said one vertical plane and the auto-correlation of said traces 
lying in said another vertical plane are obtained to normalize said cross- 
correlations in the in-line direction and in the cross-line direction. 

5 

6. The method of Claim 1 , wherein step (c) comprises the steps of 
computing the zero mean lagged cross-correlation in said in-line direction, 
and computing the zero mean lagged cross-correlation in said cross-line 
direction. 

10 

7. The method of Claim 6, wherein step (d) comprises the steps of 
identifying the most positive value of said zero mean lagged cross- 
correlation in said in-line direction, and identifying the most positive value of 
said zero mean lagged cross-correlation in the cross-line direction 

15 

8. The method of Claim 7, wherein step (d) comprises the step of 
computing the geometric mean between said two most positive values. 

9. The method of Claim 1 , where in step (a) said set of seismic 
20 signal traces comprises a plurality of amplitude-versus-horizontal 

coordinate-and-time traces of seismic data. 

10. The method of Claim 1 , where in step (a) said signal traces are 
digitally formatted. 

25 

1 1 . The method of Claim 1 , wherein each of said horizontal slices 
of step (b) extends over about 100 milliseconds. 

12. A method of locating subterranean features, faults, and 
contours, comprising the steps of: 

a) obtaining 3-D seismic data covering a pre-determined 
volume of the earth; 

b) dividing said volume into an array of relatively small 
three-dimensional cells wherein each of said cells is characterized by at 
least three laterally separated and generally vertical seismic traces located 
therein; 

c) measuring in each said cells the coherency/similarity of 
said at least three traces relative to two pre-determined directions; and 
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d) displaying said coherency/similarity of said cells to form 
a two-dimensional map of subterranean features. 

13. The method of Claim 12, where in step (c) said pre-determined 
5 directions are mutually perpendicular; and wherein said coherency/similarity 
of said cells is measured as a function of the cross-correlation between two 
traces in one direction and the cross-correlation between two traces in a 
direction that is perpendicular to said one direction. 

10 14. The method of Claim 13, wherein said coherency/similarity of 

said cells is measured as a function of the greatest cross-correlation in each 
of said two directions. 

15. The method of Claim 14, wherein said coherency/similarity is 
1 5 proportional to the geometric mean of said two greatest cross-correlations. 

16. The method of Claim 12, further including the step of: 

e) displaying the coherencies/similarities of successive 
vertically separated horizontal arrays of three-dimensional cells to identify 

20 relative space and time invarient features. 

17. The method of Claim 12, further including the steps of: 

e) selecting at least one two-dimensional pre-determined 

contour; 

25 f) dividing said contour into at least one flat 

two-dimensional surface; 

g) locating said one two-dimensional surface within a 
plurality of three-dimensional groups of said cells; 

h) measuring the median of the coherency/similarity of 
30 each group of cells relative to said one two-dimensional surface; 

i) relocating said one two-dimensional surface in said 
group of cells and remeasuring said median of said coherency/similarity to 
find a location in each group where said median is a maximum; and 

j) assigning said maximum coherency/similarity of each 
35 group of cells to one location within said group of cells. 
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18. The method of Claim 1 7, further including the steps of: 

k) obtaining the difference between the values obtained in 
step (c) and step (j): and 
5 I) displaying said difference in coherency/similarity values 

to identify subterranean shapes generally the same as said pre-determined 
contour in preference to other shapes. 

19. A method of prospecting for hydrocarbon deposits, comprising 
10 the steps of: 

a) obtaining stacked seismic signal amplitude data traces 
over at least three pre-determined three-dimensional volumes of the earth, 
at least one of said volumes having a known hydrocarbon deposit; 

b) dividing each of said volumes into an array of relatively 
15 small three-dimensional cells, wherein each of said cells has at least three 

laterally separated seismic traces located therein; 

c) measuring in each of said cells the coherency values of 
said seismic traces; 

d) comparing coherency values of said one volume having 
20 a known hydrocarbon deposit to the coherency values of the other 

three-dimensional volumes; and 

e) drilling into a location in one of said two volumes having 
coherency values comparable to the coherency values of said one volume. 

25 20 - TnG method of Claim 19, where in step (d), the most positive 

coherency values are identified and compared. 

21 . The method of Claim 1 9, wherein step (c) includes the steps of: 

0) measuring in each of said cells the cross- 
30 correlation between one pair of traces relative to one vertical plane to obtain 
an in-line value and measuring the cross-correlation between another pair 
traces relative to another vertical plane to obtain a cross-line value; and 

(ii) combining said in-line value and said cross-line 
value to obtain a coherency value for said cell. 

35 

22. The method of Claim 21, wherein said coherency value is 
assigned to the center of said cell. 
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23. The method of Claim 19, where in step (d) in each volume the 
coherency values assigned to the cells of one generally horizontal plane are 
compared to the coherency values of an adjacent horizontal plane. 

5 24. The method of Claim 19, wherein step (c) includes the steps of 

computing normalized cross-correlations of adjacent seismic traces in each 
of said cells, and identifying the greatest of said normalized cross- 
correlations. 



10 25. In seismic exploration wherein the amplitudes of reflected 

seismic energy are recorded as a function of time to produce a series of 
seismic traces, a method comprising the steps of: 

(a) acquiring 3-D seismic data over a predetermined 
volume having at least one of a geological formation, fault and an erosional 

1 5 unconformity; 

(b) comparing the similarity of nearby regions of seismic 
data of said volume by: 

(i) determining the zero mean lagged cross- 
correlation in the in-line direction between a first trace and a 

20 second trace and determining the zero mean lagged cross- 

correlation in the cross-line direction, between said first trace 
and a third trace, 

(ii) identifying the most positive of said zero mean 
lagged cross-correlations in said in-line direction and in said 

25 cross-line direction; 

(iii) determining the geometric mean of said most 
positive zero mean lagged cross-correlation in the in-line 
direction and said zero mean lagged cross-correlation in the 
cross-line direction; and 

30 (iv) repeating steps (i) through (iii) for substantially all 

of the traces comprising said 3-D seismic data; and 

(c) displaying said geometric mean values to identify at 
least one of a geological formation, a fault and an erosional unconformity. 

26. The method of Claim 25, further including the steps of: 
applying a median planar operator to a plurality of said geometric mean 
values to find the dip and dip azimuth of a plane that best aligns with a zone 
of high signal discontinuity within said volume. 
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27. A method of seismic exploration, comprising the steps 

of: 

a) obtaining a 3-D seismic data set comprising seismic 
5 signal traces distributed over a volume of the earth; 

b) dividing said volume into spaced apart horizontal slices 
and dividing each slice into cells that are arranged into laterally extending 
rows and columns, each of said cells having three seismic traces that 
generally extend vertically through said cells and in two generally vertical 

1 0 mutually perpendicular planes; 

c) computing in each cell the cross-correlation between 
said traces lying in said one vertical plane to obtain an in-line value and the 
cross-correlation between said traces lying in said another vertical plane to 
obtain a cross-line value; 

1 5 d) computing a coherency value for each of said cells, said 

coherency value being representative of the geometric mean of the most 
positive value of said cross-correlation in said in-line direction and the most 
positive value of said cross-correlation in the cross-line direction; and 

e) displaying said coherency values of said cells across at 

20 least one of said horizontal slices. 
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d) displaying said coherency/similarity of said cells to form a 
two-dimensional map of subterranean features. 



5 13. The method of Claim 12, where in step (c) said pre-determined 

directions are mutually perpendicular; and wherein said coherency/similarity 
of said cells is measured as a function of the cross-correlation between two 
traces in one direction and the cross-correlation between two traces in a 
direction that is perpendicular to said one direction. 

o 

14. The method of Claim 13, wherein said coherency/similarity of 
said cells is measured as a function of the greatest cross-correlation in each 
of said two directions. 



15 15. The method of Claim 14, wherein said coherency/similarity is 

proportional to the geometric mean of said two greatest cross-correlations. 

16. The method of Claim 12, further including the step of: 

e) displaying the coherencies/similarities of successive 
20 vertically separated horizontal arrays of three-dimensional cells to identify 
relative space and time invarient features. 

17. In seismic exploration wherein 3-D seismic data comprising 
reflected seismic energy is recorded as a function of time to produce a series of 

25 seismic traces and wherein a computer is used that is adapted to process such 
seismic traces, an article of manufacture comprising: 

a medium that is readable by a computer and that carries instructions 
for said computer to perform a process comprising the steps of: 

(a) accessing 3-D seismic data over a predetermined 

30 volume; 

(b) comparing the similarity of nearby regions of said 3-D 
seismic data of said volume by: 

(i) dividing said volume into at least one horizontal slice 
and dividing said at least one slice into a plurality of cells that are arranged 
35 into laterally extending rows and columns, each of said cells having portions 
of at least three seismic traces extending therethrough including a first trace 
and a second trace that lie in one plane and a third trace that with said first 
trace lies in another plane that is gen rally at right angles to said one plane: 
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(ii) m asuring across each of said cells the coherency of 
said traces lying in said one plane to obtain a first coherency value and 
measuring the coherency of said traces lying in said another plane to obtain a 

5 second coherency value; 

(iii) combining said first coherency value and said second 
coherency value in each cell to obtain one coherency value that is 
representative of the coherency of said seismic traces in said cell; and 

(iv) storing said one coherency value of each cell of said 
10 at least one horizontal slice in a form for displaying said combined coherency 

values as a seismic attribute map . 

18. The article of manufacture of claim 17. wherein said medium 
carries instructions for said computer to perform step (ii) by determining the 

15 cross-correlation between said first trace and said second trace, and 
determining the cross-correlation between said first trace and said third trace. 

19. The article of manufacture of claim 18, wherein said medium 
carries instructions for said computer to perform step (ii) by determining the 

20 zero mean lagged cross-correlation between said first trace and said second 
trace, and determining the zero mean lagged cross-correlation between said 
first trace and said third trace. 

20. The article of manufacture of claim 19, wherein said medium 
25 carries instructions for said computer to perform step (iii) by identifying the 

most positive of each of said zero mean lagged cross-correlations. 

21. The article of manufacture of claim 20, wherein said medium 
carries instructions for said computer to perform step (Hi) by further 

30 determining the geometric mean of said two most positive zero mean lagged 
cross-correlations. 



22. A method of seismic exploration, comprising the st ps of: 
35 a ) obtaining a 3-D s ismic data set comprising seismic 

signal traces distributed over a volume of the earth; 

b) dividing said volume into spaced apart horizontal slices 
and dividing each slice into cells that are arrang d into laterally extendinq 
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rows and columns, each of said cells having at least three seismic traces that 
generally lie in two generally vertical and mutually perpendicular planes; 

c) computing in each cell the cross-correlation between said 
traces lying in said one vertical plane to obtain an in-line value and the cross- 

5 correlation between said traces lying in said another vertical plane to obtain a 
cross-line value; 

d) computing a coherency value for each of said cells, said 
coherency value being representative of a mean of the most positive value of 
said cross-correlation in said in-line direction and the most positive value of 

10 said cross-correlation in the cross-line direction; and 

e) displaying said coherency values of said cells across at 
least one of said horizontal slices. 

23. An apparatus, comprising: 

15 a) recorded means readable by a computer and carrying 

instructions for a process comprising the steps of: 

(1) accessing data comprising seismic signal traces 
distributed over a pre-detenmined three-dimensional volume of the earth; 

(2) arranging said three-dimensional volume into a 
20 plurality of vertically stacked horizontal slices and arranging at least one of 

said slices into a plurality of cells that are arranged into laterally extending 
rows and columns, each of said cells having portions of at least three seismic 
traces located therein, each of said portions of said traces generally 
extending through said cells, and a first trace and a second trace in said cell 
25 lying in one plane and a third trace and said first trace in said cell lying in 
another plane that is generally at an angle to said one plane; 

(3) computing in each of said cells the cross- 
correlation between said traces lying in said one plane to obtain an in-line 
value and computing the cross-correlation between said traces lying in said 

30 another plane to obtain a cross-line value; and 

(4) combining said in-line value and said cross-line 
value to obtain a coherency value for each of said cells. 

24. Th apparatus of Claim 23, wherein step (3) further includes the 
35 steps of: obtaining the auto-correlation of said trac s lying in said one plane; 

and obtaining the auto-correlation of said traces lying in said another plane to 
normalize said cross-correlations in said in-line direction and in said cross-line 
direction. 

AMENDED SHEET (ARTICLE 19) 



WO 96/18915 



PCT/US95/13644 



21 

25. The apparatus of Claim 23, wherein step (4) comprises the 
steps of: computing the zero mean lagged cross-correlation in said in-line 
direction; and computing the zero mean lagged cross-correlation in said 

5 cross-fine direction. 

26. The apparatus of Claim 25, wherein step (4) comprises the 
steps of: identifying the most positive value of said zero mean lagged cross- 
correlation in said in-line direction; and identifying the most positive value of 

10 said zero mean lagged cross-correlation in said cross-line direction. 

27. The apparatus of Claim 26, wherein step (4) comprises the step 
of computing a mean of said two most positive values. 

15 28 - ,n a computer adapted to receive 3-D seismic data and having a 

display for depicting processed 3-D seismic data, an article of manufacture 
comprising: 

a) a medium that is readable by the computer and that 
carries instructions for the computer to perform a process comprising the 
20 steps of: 

(1) arranging the 3-D data into a plurality of cells that 
are arranged into laterally extending rows and columns, each of said cells 
having portions of at least three seismic traces located therein, including a 
first trace and a second trace that lie in one plane and including a third trace 

25 that with said first trace lies in another plane that is at an angle to said one 
plane; 

(2) calculating in said cells representations of the 
cross-correlation between said traces lying in said one plane and calculating 
representations of the cross-correlation between said traces lying in said 

30 another plane; and 

(3) combining said representations of the 
cross-correlation between said traces lying in said one plane and said 
representations of the cross-correlation between said traces lying in said 
another plane to obtain one coherency value for ach of said c lis. 



35 



29. The article of manufacture of Claim 28, wherein said medium 
carries instructions for the computer to perform step (3) by combining a 
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representation of the greatest cross-correlation in said one plane and a 
representation of the gr at st cross-correlation in said another plane. 

30. The article of manufacture of Claim 29, wherein said medium 
5 carries instructions for the computer to combine said representations of the 
greatest cross-correlation in said one plane and the greatest cross-correlation 
in said another plane by calculating a representation of a mean of said two 
greatest cross-correlations. 

10 31. In a computer having stored therein 3-D seismic data that 

covers a pre-determined volume of the earth, a device that is readable by the 
computer and that carries instructions to perform a process comprising the 
steps of: 

(1) digitally sorting said data into an array of relatively 
15 small three-dimensional cells wherein each of said cells is characterized by at 

least three laterally separated and generally vertical seismic traces located 
therein; 

(2) calculating in each said cells a coherency value 
from said at least three traces relative to two pre-determined directions; and 

20 (3) storing said coherency values of said cells for the 

computer to display a two-dimensional map of subterranean features 
represented by said coherency values. 

32. The device of Claim 31, where in step (2) said two pre- 
25 determined directions are mutually perpendicular; and wherein each 

coherency value is calculated as a function of the cross-correlation between 
two traces in one of said two mutually perpendicular directions and the cross- 
correlation between two traces in the other of said two mutually perpendicular 
directions. 

30 

33. The device of Claim 32, wherein said coherency value of step 
(2) is computed as a function of the greatest cross-correlation in said one 
direction and the greatest cross-correlation in said other direction. 

35 34. The device of Claim 33, wherein said coherency value of step 

(2) is a function of the geometric mean of said two greatest cross-correlations. 
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35. A method of prospecting for hydrocarbon deposits, comprising 
the st ps of: 

a) obtaining 3-D seismic data over a pre-determined three- 
dimensional volume of the earth; 
5 b) using a computer and a program for said computer that 

instructs said computer to perform the following steps: 

(1) reading said data and dividing said volume into an 
array of relatively small three-dimensional cells, wherein each of said cells 
has at least three laterally separated seismic traces located therein; and 
10 (2) calculating in each of said cells coherency values 

of said seismic traces; and 

c) using said computer to display said coherency values; 

and 

d) using said display to identify geological features and 
15 locations that are indicative of the location of an oil or gas deposit. 

36. The method of Claim 35, further including the step of 

e) drilling a well at a location identified in step (d). 

20 37. The method of Claim 35, wherein step (2) is performed by: 

(i) measuring the cross-correlation between one pair of 

traces relative to one vertical plane to obtain an in-line value and measuring 

the cross-correlation between another pair traces relative to another vertical 

plane to obtain a cross-line value; and 
25 0*0 combining said in-line value and said cross-line value to 

obtain a coherency value for said cell. 

38. The method of Claim 37, wherein step (ii) includes the steps of: 
identifying a maximum in-line cross-correlation and a maximum cross-line 

30 cross-correlation; and combining said maximum cross-correlations. 

39. The method of Claim 38, wherein said maximum cross 
correlations are combined by computing their mean. 

35 40. A device, comprising: 

a) a recording that is readable by a computer and that 
carries instructions for a process comprising th steps of: 
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(1) reading data representative of seismic signal 
traces distributed over a pre-determined three-dimensional volume of the 
earth; 

(2) sorting said signal traces by dividing said three- 
5 dimensional volume into a plurality of relatively thin cells that are arranged 

into laterally extending rows and columns, each of said cells having portions 
of at least three seismic traces located therein, including a first trace and a 
second trace that be in one plane and a third trace that lies with said first 
trace in another plane that is at an angle to said one plane; 

10 (3) measuring in each of said cells the cross- 

correlation between said traces lying in said one plane to obtain an in-line 
value and the cross-correlation between said traces lying in said another 
plane to obtain a cross-line value; and 

(4) combining said in-line value and said cross-line 

15 value to obtain one coherency value for each of said cells. 

41. The device of Claim 40, wherein the computer includes means 
for displaying said coherency values of said cells. 

20 42. The device of Claim 40, wherein step (4) comprises the steps 

of: computing the zero mean lagged cross-correlation in said in-line direction; 
and computing the zero mean lagged cross-correlation in said cross-line 
direction. 

25 43. The device of Claim 42, wherein step (4) comprises the steps 

of: identifying the most positive value of said zero mean lagged cross- 
correlation in said in-line direction; and identifying the most positive value of 
said zero mean lagged cross-correlation in the cross-line direction. 

30 44. The device of Claim 43, wherein step (4) comprises the step of 

computing a mean between said two most positive values. 

45. A method of locating subterranean features, faults, and 
contours, comprising the steps of: 
35 a) obtaining seismic data covering a pre-determined volume 

of the earth; 
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b) dividing said volume into an array of relatively small 

LeTT™' ^ ^ ^ °' SaW °* is Notarized by atC 
three te ,era.ly separated and generally vertical seismic traces located .h^ 

traces in a direction ,ha, * perpendicular ,o said one direct" an!T ° 

46. The method of Claim 45, wherein said representation of sten mi 
,s a funct,on of the greatest cross-corre.ation in each of'said £ 2££ W 
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•s a funct.on of the geometric mean of said two cross-correlation. 
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49. The method of Claim 48, wherein step (1) is performed by 
ho , (,) measuring in each cell the cross-correlation 
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value to obtain a coherency valu for said cell. 



30 



35 



BNSDOCID: <WO 961891 5A1J_> 



AMENDED SHEET (ARTICLE 19) 



WO 96/18915 



PCT/US95/13644 



26 

50. The method of Claim 49, wherein each cell contains a plurality 
of traces in each vertical plane; wherein step (i) is performed for all traces in 
each vertical plane; and wherein step (ii) comprises the steps of: identifying a 

5 maximum in-line cross-correlation and a maximum cross-line cross- 
correlation; and combining said maximum in-line and said maximum cross- 
line cross-correlations. 

51 . A seismic map prepared by a process, comprising the steps of: 
10 (1) accessing, by means of a computer, a dataset 

comprising seismic signal traces distributed over a pre-determined three- 
dimensional volume of the earth; 

(2) dividing said three-dimensional volume into a 
plurality of vertically stacked slices and dividing at least one of said slices into 

15 a plurality of cells that are arranged into laterally extending rows and columns, 
each of said cells having portions of at least three seismic traces located 
therein, each of said portions of said traces generally extending through said 
cells, and a first trace and a second trace in each cell lying in one plane and a 
third trace and said first trace in said cell lying in another plane that is 

20 generally at an angle to said one plane; 

(3) computing across each cell the cross-correlation 
between said traces lying in said one plane to obtain an in-line value and 
computing the cross-correlation between said traces lying in said another 
plane to obtain a cross-line value; 

25 (4) combining said in-line value and said cross-line 

value to obtain one coherency value for each cell; and 

(5) displaying said coherency values of said c lis 
across at least one of said slices. 

30 52. The seismic map of Claim 51, wherein prior to step (5) said 

coherency values of said cells are digitally stored in a memory; and wherein 
step (5) is performed by printing out said coherency values in the form of an 
image representative of the subsurface. 

35 53. The seismic map of Claim 51, wherein step (3) comprises the 

steps of: computing the zero mean lagged cross-correlation in said in-line 
direction; and computing the zero mean lagged cross-correlation in said 
cross-line direction. 
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54. The seismic map of Claim 53, wherein step (4) comprises the 
steps of: identifying the most positive value of said zero mean lagged cross- 
correlation in said in-line direction; identifying the most positive value of said 

5 zero mean lagged cross-correlation in the cross-line direction; and combining 
said two most positive values. 

55. The seismic map of Claim 54, where in step (4) said two most 
positive values are combined by computing their geometric mean. 

10 
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Statement under Article 19(1) 

Proper protection of the invention requires that claims should be included not only for 
the process described by the claims but also for the steps of the process when 
recorded on computer readable media and the product (a seismic map) of the 
process. Accordingly, claims have been added to cover computer media having 
instructions for performing the process of the invention and the product of the process. 
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ARBEIT 
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VORGEHEN zutreffend, nachstehender Punkt 5 
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Dieser internationale Recherchenbericht wurde von der Internationalen Recherchenbehorde erstellt und wird dem Anmelder gemaMB 
Artikel 18 ubermittelt. Eine Kopie wird dem Internationalen Buro ubermittelt. 

Dieser internationale Recherchenbericht umfaBt insgesamt _2 Blatter. 

[X| Daruber hinaus liegt ihm jeweils eine Kopie der in diesem Bericht genannten Unterlagen zum Stand der Technik bei. 



1. 



Grundlage des Berichts 

a. Hinsichtlich der Sprache ist die internationale Recherche auf der Grundlage der internationalen Anmeldung in der Sprache 
durchgefuhrt worden, in der sie eingereicht wurde, sofern unter diesem Punkt nichts anderes angegeben ist. 



□ 



2. 
3. 



Die internationale Recherche ist auf der Grundlage einer bei der Behorde eingereichten Obersetzung der internationalen 
Anmeldung (Regel 23.1 b)) durchgefuhrt worden. 

b. Hinsichtlich der in der internationalen Anmeldung offenbarten Nucleotid- und/oder Aminosauresequenz ist die internationale 
Recherche auf der Grundlage des Sequenzprotokolls durchgefuhrt worden, das 

[ j in der internationalen Anmeldung in Schriflicher Form enthalten Ist. 

f _] zusammen mit der internationalen Anmeldung in computerlesbarer Form eingereicht worden ist. 

[ | bei der Behorde nach tragi ich in schriftlicher Form eingereicht worden ist. 

| | bei der Behorde nachtraglich in computerlesbarer Form eingereicht worden ist. 

[Q] Die Erklarung, daG das nachtraglich eingereichte schriftliche Sequenzprotokoll nicht uber den Offenbarungsgehalt der 
internationalen Anmeldung im Anmeldezeitpunkt hinausgeht, wurde vorgelegt. 

[ ] Die ErklaVrung, daf3 die in computerlesbarer Form erfa3ten Informationen dem schriftlichen Sequenzprotokoll entsprechen, 
wurde vorgelegt. 

[ ] Bestimmte Anspruche haben sich als nicht recherchierbar erwiesen (siehe Feld I). 
| | Mangelnde Einheitlichkeit der Erfindung (siehe Feld II). 

Hinsichtlich der Bezeichnung der Erfindung 

|"X~l wird der vom Anmelder eingereichte Wortlaut genehmigt. 
f~~] wurde der Wortlaut von der Behtirde wie folgt festgesetzt: 



5. Hinsichtlich der Zusammenfassung 

[Yl wird der vom Anmelder eingereichte Wortlaut genehmigt. 

wurde der Wortlaut nach Regel 38.2b) in der in Feld III angegebenen Fassung von der Behorde festgesetzt. Der 
j ] Anmelder kann der Behorde innerhalb eines Monats nach dem Datum der Absendung dieses internationalen 

Recherchenberichts eine Stellungnahme vorlegen. 

6. Folgende Abbildung der Zeichnungen ist mit der Zusammenfassung zu verdffentlichen: Abb. Nr. ] 



| ] wie vom Anmelder vorgeschtagen keine der Abb. 

[X] weil der Anmelder selbst keine Abbildung vorgeschlagen hat. 
[ ] weil diese Abbildung die Erfindung besser kennzeichnet. 
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A. KLASSIFIZIERUNG DES ANMELDUNGSGEGENSTANDES 

IPK 7 G01V1/28 



Nach der Internationalen Pate ntklassifikat ion (IPK) Oder nach der nationalen Klassifikation und der IPK 



B. RECHERCHIERTE GEBIETE 



Recherchierter Mindestprufstoff (Klassifikationssystem und Klassifikationssymbole ) 

IPK 7 G01V 



Recherchierte aber nicht zum Mindestprufstoff gehorende Veroffentlichungen, soweit diese unter die recherchierten Gebiete fallen 



Wahrend der internationalen Recherche konsultierte elektronische Datenbank (Name der Datenbank und evtl. verwendete Suchbegriffe) 

EPO-Internal , INSPEC, COMPENDEX 



C. ALS WESENTLICH ANGESEHENE UNTERLAGEN 



Kategorie 0 Bezeichnung der Veroffentlichung, soweit erforderlich unter Angabe der in Betracht kommenden Teile 



Betr. Anspruch Nr. 



WO 96 18915 A (AMOCO CORP) 
20. Juni 1996 (1996-06-20) 
in der Anmeldung erwahnt 
das ganze Dokument 

US 5 706 194 A (NEFF DENNIS B 
6. Januar 1998 (1998-01-06) 
Spalte 3, Zeile 12 - Zeile 48 



ET AL) 



W0 97 13166 A (AMOCO CORP) 

10. April 1997 (1997-04-10) 

in der Anmeldung erwahnt 

Seite 4, Zeile 10 -Seite 5, Zeile 5 

WO 97 39367 A (AMOCO CORP) 
23. Oktober 1997 (1997-10-23) 
in der Anmeldung erwahnt 
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Weitere Veroffentlichungen sind der Fortsetzung von Feld C zu 
entnehmen 



Siehe Anhang Patentfamilie 



° Besondere Kategorien von angegebenen Veroffentlichungen 

"A' Veroffentlichung, die den allgemeinen Stand der Technik definiert, 
aber nicht als besonders bedeutsam anzusehen ist 

"E* atteres Dokument, das jedoch erst am oder nach dem internationalen 
Anmeldedatum veroffentlicht worden ist 

"L" Veroffentlichung, die geeignet ist, einen Prioritatsanspruch zweifelhaft er- 
scheinen zu lassen, oder durch die das Veroffentlichungsdatum erner 
anderen im Recherchenbericht genannten Veroffentlichung belegt werden 
soil oder die aus einem anderen besonderen Grund angegeben ist (wie 
ausgefuhrt) 

"O" Veroffentlichung, die sich auf eine mundltche Offenbarung, 

eine Benutzung, eine Ausstellung oder andere MaBnahmen bezieht 

"P" Veroffentlichung, die vor dem internationalen Anmeldedatum, aber nach 
dem beanspruchten Priorftatsdatum veroffentlicht worden ist 



■~P Spatere Veroffentlichung, die nach dem internationalen Anmeldedatum 
oder dem Prioritatsdatum veroffentlicht worden ist und mtt der 
Anmeldung nicht kollidiert, sondern nur zum Verstandnis des der 
Erfindung zugrundeliegenden Prinzips Oder der ihr zugrundeliegenden 
Theorie angegeben ist 

'X' Veroffentlichung von besonderer Bedeutung; die beanspruchte Erfindung 
kann allein aufgrund dieser Veroffentlichung nicht als neu oder auf 
erf inderischer Tatigkeit beruhend betrachtet werden 

■V Veroffentlichung von besonderer Bedeutung; die beanspruchte Erfindung 
kann nicht als auf erf inderischer Tatigkeit beruhend betrachtet 
werden, wenn die Veroffentlichung mit einer oder mehreren anderen 
Veroffentlichungen dieser Kategorie in Verbindung gebracht wird und 
diese Verbindung fur einen Fachmann naheliegend ist 
Veroffentlichung, die Mitglied derselben Patentfamilie ist 
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